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Abstract. In this paper, we suggest a new heterogeneous multiscale method (HMM) for 
the time-harmonic Maxwell equations in locally periodic media. The method is constructed by 
using a divergence-regularization in one of the cell problems. This allows us to introduce fine- 
scale correctors that are not subject to a cumbersome divergence-free constraint and which can 
hence easily be implemented. To analyze the method, we first revisit classical homogenization 
theory for time-harmonic Maxwell equations and derive a new homogenization result that 
makes use of the divergence-regularization in the two-scale homogenized equation. We then 
show that the HMM is equivalent to a discretization of this equation. In particular, writing 
both problems in a fully coupled two-scale formulation is the crucial starting point for a 
corresponding numerical analysis of the method. With this approach we are able to prove 
rigorous a priori error estimates in the H(curl)- and the H”^-norm and we derive reliable and 
efficient localized residual-based a posteriori error estimates. 
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1 Introduction 

The behavior and propagation of electromagnetic fields is studied in many phys¬ 
ical applications, for instance in the large area of wave optics. Periodic and 
locally periodic materials are considered with growing interest, for example in 
the application of photonic crystals (see [5^ for an introduction), as they can 
show unusual behavior, such as photonic band gaps and even negative refraction 
(see e.g. [niiHiiniiss])- However, a thorough mathematical understanding of 
these phenomena is still lacking. Therefore, one major goal is to develop effi¬ 
cient numerical schemes to simulate wave propagation in periodic materials and 
to rigorously analyze the new algorithms and the errors they introduce. 

Maxwell’s equations in a linear conductive medium, subject to Ohm’s law, 
are given by 

curlE(a;,t)-I-^(a;)i9tH(a;,t) =0, div(^(a;)H(a;, t)) = 0, 

curl H(a;, t) — e{x)dt'E{x, t) — a{x)'E{x, t) = j(a:, t), div(e(a;)E(a:, t)) = p{x, t). 

*This work was supported by the Deutsche Forschungsgemeinschaft (DFG) in the project 
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f Institut fiir Numerische und Angewandte Mathematik, Westfalische Wilhelms-Universitat 
Munster, Einsteinstr. 62, D-48149 Munster 
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Here, the electric field E and the magnetic field H are the unknowns, the current 
density j and the charge density p are the sources, and the permittivity £, the 
permeability p, and the conductivity a are material parameters. The time- 
harmonic Maxwell equations can be obtained by assuming periodicity in time, 
i.e. for 'll) = E,H,j,p we can make the ansatz tl){x,t) = Re(' 0 {a:)e®‘^*) with a 
complex-valued function il) and a temporal frequency w 7 ^ 0. Inserting this 
ansatz in the original equations yields the time-harmonic system 

curl(/r“^ curlE) -|- {iuja — a;^e)E = —iu)j, (1.1) 

H = curlE. (1.2) 

In this paper, we will focus on the curl-curl-problem (HI) on a bounded 
domain and with locally periodic coefficients. More precisely, let 11 C be a 
bounded domain with outer normal n on 911 and we seek E^ : H —)• with 

curl(^^^ curlE^) — k^E^ = f in H, (1.3) 

E 5 X n = 0 on dn. (1.4) 

We will assume that the real-valued parameter and the complex-valued 
parameter ks are locally periodic with periodicity length 6, where 6 is very 
small compared to 12. The boundary condition dn models the case where fi 
is surrounded by a so called perfect electric conductor. We refer to [221 [31], or 
m for a detailed motivation and further applications. 

Since a numerical treatment of (fOll-lfLl) requires discretizations with mesh 
sizes /i < (5 <C 1 , corresponding computations can easily exceed today’s available 
computer resources if tackled with a standard approach. In order to make the 
problem numerically solvable, so called multiscale methods can be applied. One 
class of multiscale methods that has been proved to be very efficient for scale- 
separated problems with local periodicity (or mild heterogeneities) is the family 
of Heterogeneous Multiscale Methods (HMM) introduced by E and Engquist 
[SIH- HMM approaches exploit structural invariants in the coefficients to 
solve local sample problems that allow to extract representative features and to 
approximate the (5-dependent multiscale solution with a computational complex¬ 
ity that is however independent of S. With this strategy the problem becomes 
solvable even for arbitrarily small values of S. Eirst analytical results concern¬ 
ing the approximation properties of the HMM for elliptic problems have been 
derived in [mu 111130]. In this contribution we formulate and analyze a new 
HMM for solving the curl-curl-problem (11.31) ^ (11.41) . 

Concerning wave propagation in general, the HMM and related multiscale 
methods for wave equations have been studied by Abdulle et.al. 013] and En¬ 
gquist, Runborg et.al. 0 [TS] [191 HO]- An HMM for the Helmholtz equation 
has been suggested in m- Eurthermore, some methods based on asymptotic 
expansions have been suggested for Maxwell’s equations, see e.g. [01133]. 

The new contribution of this article is the first formulation of a Hetero¬ 
geneous Multiscale Method for the time-harmonic Maxwell equations and its 
comprehensive numerical analysis in terms of a priori and a posteriori error es¬ 
timates. The error analysis can serve as a starting point for a locally adaptive 
version of the described HMM and the scheme itself might be applicable to other 
related problems after slight modifications. 
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The idea of the HMM is to adapt the (analytical) homogenization procedure 
to the numerical scheme. Therefore, we will first have a look at the homogeniza¬ 
tion of the time-harmonic Maxwell equations. Combining results by Wellander 
et.al. EH sni m] and Visintin |d9] . we derive a new two-scale equation for 
time-harmonic Maxwell’s equations. One essential step in the homogenization 
procedure is a divergence-regularization in order to incorporate a divergence-free 
constraint imposed on the corrector of the curl into the equation. This regular¬ 
ization also is an essential ingredient in the formulation of the new HMM. We 
will then adopt the view of the HMM as a direct discretization of the derived 
two-scale equation. This reformulation builds the crucial ingredient for an a pos¬ 
teriori analysis. It has been first developed in [30] and has then been adopted to 
other problems, as perforated domains [23] or advection-diffusion problems pSj . 
for instance. There have been several contributions on the numerical analysis for 
time-harmonic Maxwell’s equations (see the excellent book [53] for an overview, 
m for higher order finite elements, and EaiH] for a posteriori analysis). In the 
analysis of our HMM we need to combine these two approaches in a new way. 
Thereby we are able to obtain optimal error estimates. 

The article is organized as follows: In Section [5] we formulate the multi¬ 
scale curl-curl-problem and give some properties of the solution. The problem 
is homogenized with the tool of two-scale convergence in Section |21 The homog¬ 
enized formulation is the motivation and starting point for the formulation of 
the HMM in Section jT] Error estimates for this method are given in Section |S] 
All essential proofs are detailed in Section [51 

2 Problem setting 

For the remainder of this article, let H C be a bounded, simply connected 
domain with connected piecewise polygonal Lipschitz boundary dfl and outer 
unit normal n. Throughout this paper, we use standard notation: By IT*’^(H) 
we denote the space of functions on H with weak derivatives up to order I 
belonging to U‘{Vt) and we write := IT*’^(r2) for the scalar and H*(H) := 

[iJ*(H)]^ for the vector-valued case. Vector-valued functions are indicated by 
boldface letters and unless otherwise stated, all functions are complex-valued. 
The dot will denote a normal (real) scalar product, for a complex scalar product 
we will explicitly conjugate the second component by using u* as the conjugate 
complex of u. For any domain a; C we introduce the spaces 

H(curl,a;) := {u S L^(w;C^)| curlu S L^(a;;C^)} and 
H(div,a;) := {u S L^(w;C^)| divu G L^(a;;C)}. 

For a; = H we write H(curl) := H(curl, H) and H(div) := H(div, fl). These 
spaces are complex Hilbert spaces if endowed with the scalar products 

(u, v)H(curi) := / curl u • curl v* -|- u • v* dx, 

Jo. 

(u,v)H(div) := / divu divv*-b u • V* dec. 

Jq 

With the help of a trace theorem (see [53]), zero boundary values for functions 
in H(curl) can be defined as 

Ho(curl) := {v G H(curl)|v x n = 0}. 
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For higher regularity, we define for s € No the space 

H®(curl) := {u e H(curl) | u e curlu e H"(f7)}. 

Observe that H°(curl) = H(curl). Let denote the Fth unit vector in (i.e. 
{ei)j = Sij for 1 < i,j < 3). For the rest of the paper we write Y := [—i, to 
denote the S-dimensional unit cube and we say that a function v € is 

y-periodic if it fulfills v(y) = v(y + ei) for alH = 1, 2, 3 and almost every y G R^. 
With that we denote Y^(Y) := {w G v is F-periodic}. Analogously we 

indicate periodic function spaces by the subscript jj. For example, JI^(V) is the 
space of periodic -functions and we furthermore define for s € N 

fflo(V) :=!</)€ ffi(r)l£ </>(y) dy = o|. 

By X) we denote Bochner-Lebesgue spaces over the Banach space X and 
we use the short notation f{x,y) := f{x){y) for / G LP{V,;X). 

Using the above notation we make the following assumptions on the coeffi¬ 
cients. 

Assumption 2.1. The (scalar) coefficient fJ- ^ G C'°(U;L“(y)) is real-valued 
and K G {ft; L'^{Y; <C)) is complex-valued. Let ip denote Re(K), — Im(«;) or 
/i“^. Then there exist co,ci G R such that 0 < cq < 'ip{x,y) < ci for a.e. x and 
y, such that tp{-, j) is measurable for all 6 > 0 and 

lim / 'iplx, — ] dx= / tp{x,y)^ dydx. ( 2 . 1 ) 

\ dy JqJy 

Definition 2.2 (Weak solution). Define ks{x) := k{x, |), p.f^{x) := y.~^{x, j) 
and let Assumvtion \2.1\ be fulfilled. Let f G H(div) with divf = 0. We call 
G Ho(curl) a weak solution if, for all xp G Ho(curl), it fulfills 

/ yf ^{x) cnrl Els{x) ■ cml xp* (x) — ks'Es{x) ■ xp* (x) dx = / i{x) ■ xp* {x) dx. 
Jq Jq 

( 2 . 2 ) 

For fixed S, there is a unique solution to (12.2L which can be seen using 
the Lax-Milgram-Babuska theorem, [7] : Clearly, the right-hand side is a mem¬ 
ber of the dual space and the left-hand side gives a continuous sesquilinear 
form. Since Im k is bounded away from zero, we also get the coercivity estimate 
|i 35 (u, u)| > with a ^-independent constant. See [32] for the case of 

constant coefficients and m for the general computation. Hence, we also have 
the uniform estimate ||E5||H(curi) < C'llflU^ with C = (7(co,ci,U). 

In general, solutions to curl-curl-problems do not admit more than 
regularity and may have singularities near re-entrant corners of the domain, see 
m- However, if LI is convex and if /x ^,k G lU^’°°(r2 x Y), i.e. the coefficients 
are globally Lipschitz, the weak solution to ()2.2I1 fulfills E^ G H^(curl), see m- 

3 Homogenization 

As the periodicity length d is assumed to be very small in comparison to H, one 
can reduce the complexity of the problem by considering the limit d —>■ 0. This 
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process is called homogenization and can be performed with the tool of two-scale 
convergence [1] . Since the two-scale equation and the homogenized equation are 
the starting point for the construction and analysis of the numerical multiscale 
method later on, we present the essential steps in this section. 

3.1 Two-scale convergence 

Two-scale convergence is defined as (cf. |4l [27]1: 

Definition 3.1 (Two-scale convergence). A sequence (us)s>o C two- 

scale converges to a function uq G L^(fl x Y) (short form: ug uq) if it fulfills 

lim / usix)tlj (x,^] dx= uo{x,y)ii{x,y) dydx 'iiliGLf{VL-,C'((Y)). 

^ oy JqJy “ 

For more information on two-scale convergence, for example the definition 
of strong two-scale convergence, and compactness in and we refer the 
reader to I1I27], and the lecture script |23) . 

For time-harmonic Maxwell’s equations we need a two-scale convergence 
result for bounded sequences in H(curl). As H(curl) is not compactly embedded 
in (in contrast to H^), the two-scale limit in will not coincide with the 
weak limit thus making additional considerations necessary (cf. [391 p. 135]). 
We present two possible approaches, which we will combine in our analysis later 
on. With the help of the compactness theorem and integration by parts, the 
following characterization of two-scale limits can be derived (see [TTl l40l 1^ 1: 

Proposition 3.2. Let (u5)5>o C H(curl) be a bounded sequence. Then there 
exists a subsequence and functions Uq € L^(n x Y), Ui G Hjj(curl, F)) 

and (j) G q(Y)) such that 

1. U5 ^ Uo with Uo(a:) = u(a:) -|- Vy4>{x, y) and u = fy Uo(-, y) dy G H(curl), 

2. curing curl^, Uq -I-curlj^ Ui and curlcurlu in 

Using a technique called periodic unfolding, one can obtain a characteriza¬ 
tion for the curl which resembles the one in the compactness theorem (see 

my- 

2 

Proposition 3.3. Let (u5)5>o C H(curl) be a bounded sequence such that ^ 
Uq in L^(Ll X Y). Then there exists a subsequence and Ui G L^{Lt-, Hj q(F)) with 
divy Ui = 0 a.e. in LI x Y such that 

curing ^ curlu -|- curly Ui with u := J Uo{-,y) dy. 

Furthermore, it holds Uq G L^(n; Hjj(curl, F)) and we have for the whole se- 
quence 6 curing ^ curly Uq = 0 . 

We will combine this theorem with the first point of Theorem 13.21 Note 
that the condition divy Ui = 0 can be seen as a kind of gauging condition. It 
will be important for the homogenization of our curl-curl-problem, namely this 
condition will lead to the uniqueness of the two-scale solution. 
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3.2 Homogenization for time-harmonic Maxwell’s equa¬ 
tions 

In this section we present new homogenization results for the time-harmonic 
Maxwell equations in a two-scale formulation, a formulation with cell problems 
and macroscopic equations, and a corrector result. We emphasize that although 
Maxwell’s equations and curl-curl-problems have been homogenized in the liter¬ 
ature ([ 31101111 ] only to name a few), the focus has always been on macroscopic 
(homogenized) problems as (13.71) . but not on two-scale limit equations. 

Theorem 3.4 (Two-scale equation). Under the same assumptions as in Def- 
inition fOl let G Ho(curl) be a solution of (12.211 . Then there exists a so¬ 
lution triple (E, Ki, K 2 ) of functions E G Ho(curl), Ki G L^(r2; Hj q(F)) with 
divy Ki = 0 a.e., and K 2 G Hj q{Y)) such that 

E 5 ^ E OT Ho(curl), EiE-I-Vyi 4 r 2 , curl E 5 curl E-|-curIj, Ki. 

Considered m Ho(curl) xL2(q; Hj_^(r)) the triple (E,Ki,i4:2) 

is the unique solution of 


IJ^ ^{x,y){cuT\E{x)cuT\yKi{x,y)) ■ {cut\xI}*{x) curly 

Ja Jy 

divy Ki(a;, j/) div^ xf>\{x,y) 

- K.{x,y){E{x) + \IyK 2 {x,y)) • (V’*(a:) -H Vy-if^ix^y)) dydx 
= f f(a:) ■ ip*{x) dx 

Jn 


VV’ G Ho(curl),V’i G L^{n;Ul,{Y)),rb2 G L\n;Hl,{Y)). 


(3.1) 


The proof is postponed to Section 13 but let us name the three important 
steps it consists of. First, using the two-scale convergence results mentioned 
above and inserting a special test function, we obtain a two-scale equation sim¬ 
ilar to (ED, but with divergence-free constraint. Second, we incorporate the 
divergence-free constraint into the equation (see also the remarks below). Third, 
the uniqueness of the two-scale solution is shown, which gives the convergence 
for the whole sequence of solutions E 5 . 

Remark 3.5 (Divergence-regularization). In order to determine Ki in the two- 
scale equation, one has to find u G H(curl, F) n H(div, F) with divu = 0 a.e. 
in Y and such that for all xf G H(curl, F) 


fi ^ curl u • curl if* dy = 0 , 


(3.2) 


and with appropriate boundary conditions, here periodic ones (for Ki, there is 
also a right-hand side, which we do not consider for simplicity). The divergence- 
free constraint div u = 0 is necessary to guarantee the uniqueness of a solution, 
as otherwise the solution is only determined up to a gradient term. However, the 
divergence-free constraint causes some problems in the implementation of corre¬ 
sponding numerical methods, as e.g. divergence-free finite elements are difficult 
to find or to implement. 
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With divergence-regularization, we now look for u G H^(y) such that 


JY 


Clearly, any divergence-free solution to (1^ will also solve (ESI). On the other 
hand, we can insert a test function ■0 = 'V(j) for (j) G ff^(Y) and obtain 


j div u A(j)* dy = 0. 


If the domain is convex, there is 4> G H'^{Y) with A(p = divu due to elliptic 
regularity theory, since divu G LfifY). Hence, divu = 0 almost everywhere. 

The convexity of the domain is an essential assumption in the divergence- 
regularization, so that the method cannot be applied on arbitrary domains. We 
emphasize that divergence-regularization is only needed in the corrector here and 
the corresponding problems are always posed on the unit cube Y or at most on a 
parallelepiped. Thus, convexity is guaranteed and does not impose any additional 
constraint. 

There are other possibilities to deal with a divergence-free constraint. The 
introduction of Lagrange multipliers (see \1^ ) leads to a mixed problem, which 
increases the computational costs and complicates the error analysis. The s- 
regularization suggested in m makes the reformulation of the HMM later on 
(Proposition \4.f^ impossible, since different orders of derivatives appear. Thus, 
we choose divergence-regularization as it easily gives an eguivalent problem, pre¬ 
serves coercivity, and can be implemented in the HMM framework as well. 

Definition 3.6 (Cell problems and homogenized matrices). The cell problems 
are to find functions G L^(n; Hj Q(y)), Vk G Hj q{Y)) so that a.e. in 

LI there holds 


1 yL~'^{x,y){ek -f curly Vfe(x, y)) • curl0*(y) 


(3.4) 

-1- divy Vfc (x, y) div 0* (y) dy = 0 

V0 G Hj,o(y), 


1 K{x,y){ek -f VyVk{x,y)) ■ V0*(?/) dy = 0 

V0 G Hl^{Y). 

(3.5) 


With the (unique) solutions of the cell problems (|3.4II - (13.51) we define the ho¬ 
mogenized matrices 




M ^{x,y){S^k + {cm\yVk{x,y))i)dy, 


^hom 

^i,k 



K{x,y){Sik + {\7yVk{x,y))i)dy, 


i,k = 1, 2, 3. 


(3.6) 


The homogenized matrices are used to formulate the macro-scale problem 
for E. It has the same structure as our original problem except that the material 
parameters are now matrices and no scalar functions. 


Theorem 3.7 (Equivalence of two-scale and homogenized equation). The triple 
(E,Ki,i4r2) is the unique solution of (13.111 iffEG Ho(curl) solves 


f 

Jn 


'curlE-curl0*-K''°'"E-0*da;= / f • 0* dx 


(3.7) 
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for all 'll) G Ho(curl) with the matrices ^ i^hom through (13.61) . 

and with correctors 'K.i,K2 defined as K.i{x,y) = X]fe=i(curlE(a:))feVfe(a:,?/), 
K2{x,y) = '^i^^i^k(x)vk(x,y), where VkjVk are solutions of the cell problems 

(Isa, disi). 

Proof. Inserting the cell problems and the definition of the homogenized matri¬ 
ces into (13.71) leads to the two-scale equation. □ 

We end this section by a corrector-type result, which relates the two-scale 
solution to the asymptotic expansion. 

Theorem 3.8 (Strong convergence in H(curl)). Let k, Ki, curia, Ki, 
curly Ki, VxK 2 and VyK2 be admissible test functions for two-scale conver¬ 
gence. Then it holds 

|E5(a;) - (E(a;) -f ,5 (Ki ( x , |) + VK2 (x, |))) 

Proof. Inserting the term in the norm into the heterogeneous sesquilinear form 
Bs, using the chain rule and two-scale convergence gives the claim. □ 

The theorem shows that the correctors Ki and K2 represent a Helmholtz 
decomposition of the first order term in the asymptotic expansion. Since on 
the gradient subspace, the H(curl)-norm and the L^-norm are equivalent, we 
see that in particular K2 carries important information about the solution E5 
itself. Thus, in contrast to the elliptic case, the correctors Ki, K2 have to be 
considered as well (and not only the weak limit E) in order to get a good L^- 
approximation of the heterogeneous solution E^. This is a crucial observation. 
Consequently, the HMM is not only constructed to approximate E, but requires 
to approximate Ki and K2 as well. 


(5->0 
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4 The Heterogeneous Multiscale Method 

The basic idea of the HMM is to use a macroscopic sesquilinear form similar to 
(1X71) for the finite element method. Instead of solving the cell problems once 
on the unit cube, local variants are set up and solved around the centers of 
the tetrahedra of some macroscopic computational grid. In order to define the 
method in more detail, let us introduce some definitions. 

Denote hy Th = {Tj\j G J} and Th = G /} conforming, shape regular, 
simplicial partitions of H and Y, respectively, where Th is additionally periodic. 
The d-scaled and Xy-shifted unit cubes are denoted by = SY -\- Xj, together 
with the mappings yj-.Yj^Y and Xj = {yj)~^:Y -^Y^. A triangulation of 
the shifted unit cubes is then given by ThfYj) = {S\S = x^j{S),S G Th}. The 
set of interior faces is defined as 

= {(j,0 G J X J\Fyi := TyDTi ^ 0,dim(F,i) = 2,j < 1} 

and £{Th) with the faces Fik is defined analogously. The direct neighbors of 
a face Fji are wf,, := Tj U T;. The neighborhoods of vertices V, faces F, and 
elements T are defined as 

UJv ■= [J Tj, UJPji ■= [J UJv, UJTj ■= [J OJy', 

jeJ,veTj v&Fji vTFj 
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and the neighborhood of neighbors of an element is defined as 


■= U U 

VGTj V'e<^v 


We define the local mesh sizes Hj := diam(Tj), hi := diam(S'i), Hji := 
diam(fj;), hik ■= diam(i<ifc), and the global mesh sizes H := maxjgjiJj and 
h := maxi^i hi. Finally, the discrete function spaces V]jq C Ho(curl) and 
C H^q{Yj) are defined as 


■= {uff e Ho(curl)|u;,|T € No VT e Th}. 

Wl{Y^) := {UH G HIo{Y^)\uh\s G P' G r^(F/)}, 


where P^ are the polynomials of maximal degree 1 and No is the lowest order 
Nedelec element of the first family, given by No := {a x cc + b|a, b G C^}. 
As in the analytical case, bold face letters indicate vector-valued functions and 
function spaces, for instance := 

With these preliminaries we can now define the HMM (see also [TOl [T71 [50] 1: 


Definition 4.1 (HMM). The HMM-approximation of (12.21) is a discrete solu¬ 
tion triple (Eff, Ri(Eff), R2(Eff)), where Eh G V|^q is defined as the solution 

of 

Bh{Eh, = if, ^h) e (4.1) 

where the discrete sesquilinear form is given by 


jeJ 


^ J ^ih ^)i{x)cui\'Ri{uH){x) ■ cmlifnix) 

- /v^(a;)R2(u//)(x) • R2('j/’ff)*(a;) dx 


(4.2) 


with the piecewise constant approximations := for 

all Si G ThfYj) and (/x“^)^ defined analogously. The local reconstructions 
Rl(Uff) G Unlyf + R2(Uff) = UHixk)\Yf + ^ yUh with Uh G WUX^) 

are defined as the solutions of the local cell problems 


/ {fj, ^)^(a:)curlRi(ui/) •curlt/;);-bdiv(Ri(uff) 

Jy^ 

3 

ViA, G ^liY^^)), 

f ni{x)R2{uH)-Yrhdx = 0 Vi^h€W^iYf). 

JY^ 


Uh) div i/jh dx = 0 


We now reformulate the reconstructions of the HMM solution triple to draw 
a parallel between them and the analytical correctors. 

Remark 4.2 (Role of the reconstructions). Let (Ejj, Ri(Eij), R2(E^f)) denote 
the HMM-approximation from Definition \4-. 1\ Setting Kj^i(E//) = R-j.i(Eff) - 
Eh, we have Ejp G Furthermore, denote by Kp 2 {Eh) G W^{Yj) 

the function fulfilling NKj^ 2 iEH) = Rj,2(Eij) — EHixj). We then define the 
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discrete fine-scale corrections S S'ij(fi; W^(y)) and Kfi, 2 (EiH) S 

SH{n;Wl{Y)) as 

K,,i(Eff)(a;,y)|T,xY := i(E^)(5y), 

Kh,2(E,H){x,y)\TjxY ■= ■^^j,2(E//)(5y), 
where the space of x-constant discrete functions is defined as 

Sh[VL-WI{Y)) := {uh e L^[VL,Hl^{Y))\uh{;y)\T, e P°Vj G J,?/ € F 
and Uhix, •) G Wl^{Y) Vx G ^i}. 

T/ie discrete fine-scale corrections K;i4(Eij), Kfi, 2 (EH) are discrete counter¬ 
parts of the analytical correctors Ki and introduced in Theorem \3.4\ The 
specific relation of both will be clear from Proposition \4.3\ below. Therefore, 
these corrections (or equivalently the reconstructions) are an important part of 
the HMM-approximation. As discussed at the end of Section the correctors 
carry important information on the solution and cannot be neglected as higher 
order terms (in contrast to the elliptic case). In form of the fine-scale correc¬ 
tions, the observation transfers to the numerical scheme and the discrete setting. 

Having observed this correspondence, we can now reformulate the whole 
HMM to see that it is a direct discretization with numerical quadrature of the 
two-scale equation iH). See [30] for the approach in the elliptic case. 

Proposition 4.3 (Reformulation of the HMM). Define the piecewise constant 
approximations kh on D x Y by Kh{x,y)\TjxSi ■= K{xj,yi) and in the 
same way. Furthermore, let ^h,i, Kh ,2 be the discrete fine-scale corrections as 
definedThen (E^^, K^,,i(E^^),i^;,,2(Eff)) G x Sh{^-,^1{Y)) x 

Sh{^',W(,{Y)) is a solution of 


h'h^{x,y){cur\EH{x) 


■ curlj^ 


L (Eff )(x, y)) ■ {'ipnix) curly 'ipl{x, y)) 


-f divy Kh,i{'EH){x, y) divy ipHx, y) 

- k{x, ?/)(Eh(x) -f YyKh,2i'EH){x, y)) • (^*{x) -f Yytpl{x, y)) dydx 



G X L^{n;Wi{Y)) x L^in-,W(,iY)). 


Proof. We treat the two terms of the sesquilinear form separately, but with 
basically the same procedure. For the first term we see from the definition of 
the reconstruction Ri that for all ■0^ G '^l(Yf)) it holds 


0 = 


(/i ^)^(x) curfo(E//-f Rj,i(Eh) 


Eh)(x) ■ curR0^(x) 


+ div2;(Ry,i(E/f) - Eh)(x) diVa; 'ifl(x) dx. 

Using the transformation formula and writing x = Xj{y), we derive 

0 = S^J p~^(xj, ^^^^curR0^(x^(j/))-(curfo E//(xy)-hcurR Ky,i(Eif)(Xy (y))) 
-f divj; Kj,i(Eif)(x^(y)) div^, xfl^x^iy)) dy, 
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since curl^; E//(a::) is constant on each Tj. Using the definition of Kft^i(E//) and 
defining G W^(y) as tl’kiu) = S^t with the chain rule 

4(y)\ , 


Jy^ -^) curly 

• (^curU Eif(xj) + curly K/,,i(E/f) 


/ x^Av)\ ~* 

+ divy divyi/j^(^-^j dy. 

As the integrand is U-periodic and — = y+^, we finally obtain the Galerkin 
orthogonality 

0 = J yJ^^{xj,y)icuTlxEHixj) + curly Khx{'E^H){xj,y)) • curly ^^(y) 

5 (! 

+ divy Kh,i(EH)ixj,y) divy ■0^(y) dy. 

In the same way, we deduce 

y curU IIi{Eh){x) ■cur\^pH{x) dx 


Xjiy) 

Hy)' 


X 

* /X 


3^J 


'^\T3\ j y curiv^^(a;y) 

teJ 

• (curlE_f/(a;y) + curly 'Kh,i[EH){xj,y)) dy 

■ / / MftHa;,2/)(curlE//(a;) + curly K;i,i(EH)(a;,y))-(^^(a;) + curly V’ii;(a:,y)) 
JQ, Jy 

+ diVy Kh,i(EH)ix, y) diVy -iplix, y) dydx. 


In the last equality we used that the given quadrature rule is exact for the 
integrands and we employed the Galerkin orthogonality. For the second term 
in the sesquilinear form, one can perform the same steps to reformulate the 
problem. For this term the computations are very similar to the elliptic case 
discussed in [3ni Lemma 3.5]. □ 


Conclusion 4.4. Let us note that the result of Theorem \S.8\ is still valid if we 
replace VK 2 by S~^'VyK 2 . This implies that we can approximate E5 in H(curl) 
by E( x) + 5Ki (x, I) + VyK 2 (x, fj. Consequently, exploiting Proposition EH 
we see that our final HMM-approximation Ehmm to E^ can be written as 

Ehmm(3^) ^h{^) + SKhpiEn) [x, I) + yyKH, 2 iEH) [x, . 


5 A priori and a posteriori error analysis 

Based on the reformulation of the HMM in Proposition 14.31 we will give the 
main a priori and a posteriori error estimates in Theorems l5.2l[S75115.51 and 15.71 
All error estimates will be derived in the ” two-scale energy norm” 

|1(U, Ui,U2)||e(GxiJ) 

:= II curlu -I- curly Ui||i2(GxiJ) + II divy Ui \\l‘^(gxR) + I|u -f VyU2||L2(Gxfi:) 


II 









for G X i? C fl X Y an open subdomain. If the norm is to be taken over 
n xY, we will just write || • ||e. Let us furthermore define the error terms eo := 
E —E//, d := Ki —K/j 4 (E//), and 62 := K 2 —Kh, 2 {^H)- We will only estimate 
these errors and leave the modeling error E 5 — (E + 5 (Ki (•, + VLf 2 (•, 7 )))> 

introduced by homogenization, apart (cf. Theorem 13.81 and Conclusion 14.411 . 

Assumption 5.1. On top of the periodicity of the coefficients, we also assume 

X Y), 

i.e. the coefficient functions are globally Lipschitz, and 0, is a convex domain. 
This assumption will be required for the a priori estimates f Theorems \ , 5.31 and 
15.5j) . but not for the a posteriori estimates fTheorems \ 5. 5\ and \5. Tp . 

Theorem 5.2 (A priori estimate in the energy norm). Under Assumvtions 
and fXH the following a priori estimate for the error between the homogenized 
solution and the HMM-approximation resp. their correctors holds: 

||(eo,ei,e 2 )|U < C{H + /i)||f|| l 2 (q) with C = 

Theorem 5.3 (A priori error estimate with dual problems). Under the same as¬ 
sumptions as in Theorem \5.3l the Helmholtz decomposition of the error between 
the continuous solution E and the HMM-approximation Yih 

E - E// = V 6 » + z with 9 e Hi (fl), zAViLg 

satisfies 

ll^llL2(n) + ||z||L2(f2) < C{H^ + /l^)||f||i2(f2) + Ct] approx l|f||L2(n), 

where rjapprox = \\k-K h\\L'=-(axY)} is a data approx¬ 

imation error arising from numerical quadrature. The constant C only depends 
on the domain 12, the coefficients and n, but not on the periodicity param¬ 
eter 5 or the mesh sizes. 

Remark 5.4. In the elliptic case, the L^-norm of the error converges with 
quadratic rate. This better convergence is obtained by posing a dual problem 
and using the Aubin-Nitsche trick. The above theorem shows how the result can 
be transferred to problems in H(curl); On the gradient subspace, the L'^-norm 
is of the same order as the H(curl)-norm, so that only on the complement a 
better convergence is obtained (see Remark after Thm. 49, p. 45])- Hence, 
the quadratic convergence here is (only) obtained in H~^. 

Theorem 5.5 (A posteriori error estimate). Lettn be any piecewise polynomial 
approximation oft. Under the Assumvtion \2. 1\ the error fulfills the following a 
posteriori error estimate 

||(eo, ei, 62) lU < G(^^r7|i-k 771,2) Vji.i + Vji,2) 

j&J U,i)eSiTH) 

+ C'(X! X! + vlik, 2 ') 

j&J {i,k)^e{Th) 

jeJ iGJ ie/ 
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where the constants do not depend on the mesh sizes and the periodicity param¬ 
eter 6. The local indicators are defined as 


fH+ j + yyKh,2{-,y))dy 

V],2'=Hj diva;(^J Kh{-,y){'EH + S/yKh,2{-,y))dy^ 


ul /2 

Vji,i ■■= 

ul /2 

Vji,2 ■■= H-l 


lY 

-1 






^JY 


fj,^^{-,y){cuTl'EH + cm\yKh,ii-,y)) x ndy 


Jf„ 


LHFii) 


Khi-,y)(E,H + '^yKh, 2 {-,y)) ■ ridy 


'-JY 


LHF^i 


Vj,ik,i h^i^ \\[Th (curlE// + curly K/j^i) x n + divy 
Vj,ik,2 ■= h^l ||[k/i(E// +'^yKh,2) ■ Il]_Fij.llL2(TjXFifc)’ 

Cj ■= Hj\\iH - f\\L^{Tj), 

Cji '■= WiTh^ -/r“^)(curlEij + curly K/,4)||i2(-r^.xSi) 
+ \\{F.h — k)(E,H + ^yKh,2)\\L'^{TjXSi)- 


Here and in the following, [-Jf denotes the jump across the face F. 

Remark 5.6 (Discussion of the error indicators). The error indicators can be 
split into two groups: Q and Qi are data approximation errors, which come 
from the use of numerical quadrature. The error indicators denoted by rj are 
different contributions to the discretization error: rjj^i is the element residual 
on the macro-scale, rjji are the jump residuals on the macro-scale (in normal 
and tangential direction), and rjj^ik ctre the jump residuals on the micro-scale. 
rij ^2 indicates how well the deduced equation div 2 :(KEo) = 0 is fulfilled in the 
discrete case. Here, the assumption div f = 0 (in the weak sense) has an effect 
on the error estimator: If we just assume f £ H(div), the deduced equation is 
diva;(KEo + f) = 0 and thus we have additional terms div^; in in r]j ^2 o,nd in ■ n 
^ 31,2 for the polynomial approximation of f. Furthermore, in the data 
approximation error Q we then have to take the tl{div,Tj)-norm. 

Theorem 5.7 (Lower bound on the error). With the same notations and under 
the same assumptions as in Theorem \5.5\. the following local bounds on the error 
hold: 




Vo,2 


VjiA 


Vjl,2 




1 /2 

< c{\\{eQ,ei,e2)\\e{TjxY) + Cj + )) 

i 

1 y/2 

< <^( 11 ( 60 , 61 , 62 ) ||e(T, xF) + )’ 

i 

1 /2 

< C'^IKeo, 61, e2)||e(iiF^, xY) + iCj + Cji + Ch) ) 1 

i 

1 y/2 

< C'^||(6o, 61,62)||e(cDF^, xY) + (Ea + cs) ). 

i 

< C'(||( 6 o, 6 i, 62 )||e(TiXcDp.^) + iCji + V = 1,2. 
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Thus, we have the global estimate 


iX^^li+^12)' + { Y. ^lu+4.2)^ 

3^J U,i)ee(TH) 

+ {Y Y ' 

jdJ (i,k)^e{Th) 

< c(^||(eo,ei,e2)||e + + (YY'^t) )' 

jeJ j&J lei 

Remark 5.8. Theorems \5.5\ and \5. 7| together show that the local error indicators 
are reliable and efficient with respect to the two-scale homogenized problem from 
Theorem \3.4\ Up to data approximation errors, the error and the indicators 
converge with the same rate. Thus, the indicators can be used for adaptive 
algorithms, e.g. for mesh refinement, both on the coarse and fine scale (cf. 
for related adaptive algorithms in the elliptic case). 


6 Proofs of the main results 


In this section, the essential proofs of the homogenization result (namely the 
two-scale equation) and the error estimates for the HMM will be given. 

Proof of Theorem [M As solutions to ( 1 ^ are uniformly bounded in H(curl), 
we have by Theorem 13.31 that, along a subsequence, Eq in x E))^ 

and curlEi ^ curlE -|- curly Ki. Furthermore, we have from Theorem Id.21 that 
Eq = E -I- VyK 2 , where E = fy Eo(-, g) dy is the weak limit of E^ in H(curl). 
We insert tfix) = w(a;) -b i5wi(x, f) with arbitrary w € € 

(C^(V))^) as a test function in (12.21) . Because of the assumption (12.11) 
on the parameters we can apply two-scale convergence to each of the terms in 
(1^ and thereby obtain 


Y 


fj, ^{x, y)(curl E(a;) -I- curly Ki(a;, y)) ■ (curlw*(a;) -I- curly w*(a:, y)) 
- k{x, y){E{x) -f WyK 2 {x, y)) ■ yv*{x) dydx 
f(x) • w*(x) dx. 


JQ, 


( 6 . 1 ) 


By density this also holds for test functions in Ho(curl) x L^(n; Hj q(E)). Fur¬ 
thermore, we can deduce from (EH with div f = 0 that it holds 


5 f ks{x)Es{x) ■ \/iIj*{x) dx = 0 G Hq{V.). 

Jn 


Choosing 4>{x, f) € C^(fl; C^(V)), we obtain with two-scale convergence 


/ / K(x,y)(E(x) -l-VyK 2 (x,y)) ■Vyfi*(x,y)dydx = 0. 
Jn Jy 
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Inserting this into (EU), we get (ED except for the divergence term, but with 
the additional constraint divy Ki = 0 . 

As discussed in Remark EH we can apply divergence-regularization in this 
case to obtain an equivalent problem. Looking at the method, we directly see 
that (ED is simply the regularization of ED- The equivalence of the problems 
can be seen as discussed in the remark, just insert -i/; = 0 , -02 = 0 , and 
with ip G H^q{Y)) as test function in (16.11) . 

So far we have shown ED just for a subsequence. If we can prove that 
the solution of the two-scale equation is unique, the result holds for the whole 
sequence. In fact we will prove that ED is of the form = F(^/>) with 

a continuous and coercive sesquilinear form B and a functional F. As for the 
weak solutions to (ED, Lax-Milgram-Babuska then yields the uniqueness of the 
solution. Moreover, this reformulation will be important for the error estimates 
later on. We consider the Hilbert space 

H := Ho(curl) x x L^{n;Hl,iY)) 

with its natural norm 


II(u,Ui,U2)||h — l|u||H(curl) + II *^1II L2(0;H1 (V)) + II ^2 || i (V)) ' 

Clearly, the right-hand side is in the dual space of H and the left-hand side 
defines a continuous sesquilinear form B. With the same computations as for 
the existence of a weak solution, one can also show that B is coercive with 
respect to the energy norm 


II curlu -I- curlj, Ui|||2(qxv) + II ’^illi2(nxv) + 11*^ + ^y“2||i2(nxy)- 
It remains to show the equivalence of the energy and the natural norm. It holds 
llcurlu-bcurlj^ Ui||^2(oxv) 

= / / I curlup-f I curly Uip-I-2 Re(curly Ui(a:, j/) • curlu*(a;)) dj/dx 

Jn JY 

= II curlu||^2(f2) + II curly Ui||^2(qxV) 


- 2 Re 


/ / (ui(a:, j/) X n) • u*(a:) dcrdx I 

Jn JdY ) 


=0. periodicity of ui 


= l|curlu|||2(n) + ||curlui|||2(oxF)- 

With integration by parts and the Poincare inequality we see that || curly Ui||^ 2 -|- 
II divy Ui ||^2 is equivalent to the full id^-norm. Similarly, we derive 

||U -I- VyM2|li2(f2xV) 

= l|u||i 2 (o) -b ||VyU 2 ||i 2 (oxv) + 2 Re ■nu;(x,y)dadx 

" -V- 

=0, periodicity of U2 


and again with the Poincare inequality the id^-seminorm is equivalent to the full 
dd^-norm. Due to the uniqueness of the two-scale solution the whole sequence 
converges as asserted in the theorem. □ 
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Having identified the variational formulations in Theorem 13.41 and Proposi¬ 
tion we collect some useful notation for the proofs of the error estimates in 
the next paragraphs. 

Notation 6.1 (Sesquilinear forms B, Bh and residual). We define the con¬ 
tinuous and discrete sesquilinear forms B,Bh ■ [Ho(curl) x L^(H;Hjg(y)) x 
as 


/ / /x“^(a;,y)(curlE(a;)-I-curly Ki(a;,j/)) • (curl-I-curly 

Jo, Jy 


divy Ki (x, y) divy -0*(x, y) 

- K{x,y)(E{x) VyK 2 {x,y)) ■ (0*(x) -h Vy'if 2 {x,y)) dydx, 
BhiiuH,Uh,Uh),{'lpH>'^h^’>Ph)) ■ = 


h-h^{x,y){cmluHix) -I-curly Uh{x,y)) ■ (curl0^(x) -I-curly 0))(x,y)) 


-f divy u,,(x, y) divy tplix, y) 

- i^h{x, y){uH{x) -f WyUh{x, y)) ■ {'ipnix) -t- Vyipl^x, y)) dydx. 


In addition we define the residual 

Resh : Ho(curl,H) x L"(H; Hj_o(y)) x L\n-Hlo{Y)) 

^ Ho(curl)' X L^{n; (H„-i(y))3) x L^{n; Hl^{Y)) 

as (Res,^(u,ui,u2),(0,0i,02)) := Ui, uz), ( 0 , 0 i, 02)) - (f, 0 ). 

( 6 . 2 ) 


In the following, we will write instead of Kft^i(E//), and ifh ,2 instead 
of Kh^ 2 (EiH)- If it is clear on which variables (x, y) functions depend, we will 
omit those variables for the sake of readability. C denotes a generic constant, 
independent of the mesh sizes and 5. 


6.1 Proofs of the a priori estimates 

The a priori estimates are based on the Cea lemma, dual problems and interpo¬ 
lation operators. As the assumptions on the coefficients and the domain imply 
higher (namely H^) regularity of the solution, we can use the nodal interpolation 
operators. 

Lemma 6.2 (Lagrange interpolation operator). Denote by : (7°(T) — >■ 
W^{Y) the standard Lagrange interpolation operator. Define now for v G 
L^{n;C°(Y)) by 

Ihiv){x, y) ■■= {Ihiv){x, ■)){y) - ^(/^(u)(x, •))(s) ds. 

For all V G L'^{fl;C^{Y)) fl q(Y)) this interpolation operator is well- 

defined with S Lfi{D.] W^{Y)) and the estimate 

\\Ihi.x) - u||l2(QxF) + h\itix) - v\L2(Q,^m(Y)) < Ch'^\v\L2(Q,^H'^(Y))- 

For vector functions the interpolation operator is defined component-wise and 
fulfills the same estimates. 
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□ 


Proof. For a proof we refer the reader to [53] . 

Lemma 6.3 (Edge interpolation operator). Denote by : H^(curl) ^ 
the nodal interpolation operator for the Nedelec elements. It fulfills the estimate 

I|U - df (u)llH(curl.Tj) < C'iJj||u||Hi(curl.Tj)- 

Proof. For the definition of and the proof we refer to Eg. □ 

Proof of Theorem \5.‘A Denote by 'Em G E/i G W)j(F)), and Eh G 

L^(D; VF^(y)) the unique solution of 

B{(EH,Eh, Eh),{^lpH,'4>h:'‘Ph)) = 

G G L\n;Wl{Y)),fh G L^n-WliY)). 

According to Cea’s Lemma it holds 
||(E-E^^,Ki-E^,iF2-E^)||e 

< C( inf ||E — ■0^||H(curl) + illf_ |Kl — 

V.;.eL2(n;Wi(Y)) 

+ inf_ \K2 - tph\L^(n-,{¥)))■ 

V-hGL2(n;lYi(>")) 

With the interpolation estimates of Lemmas 16.21 and 16.31 we now derive 

||(E-Eff,Ki-E;„iF2-Eft)||e 

< C(||E - /i(E)||H(curi) + |Ki - 4 ^(Ki)U2(0;HMY)) 

+ \K2 - ltiK2)\L2(n.m{Y))) 

< C'(i?||E||H(curl) + ^|Ki|L2(n;H2(Y)) + h\K2\L^(n.H^(Y))) 

<CiH + h)\\{\\L^n), 

where in the last inequality we used regularity and stability results for the 
analytic solution. (Note that because of our assumptions on the parameters 
and on fl the two-scale solution admits regularity.) Furthermore, because 
of the definition of (E/f,E^,i?/i) it holds 

\\(Eh - Eh, Eh - K;,,i(Eh), Eh - iF,,.2(E^^))||2 

< C\{Bh - B){{EH,KhA^H),Kh,2{EH)), 

{Eh - Eh, Eh — Kft4(E//), Eh — Kh,2{EH)))\ 

< Cm.ayi{\\p.f^ - M”^llL“(nxY), ll^/i - K||L“(axY)} 

||(E_f/, K;i4(E/f), Kh^2{EH))\\e 

||(Ejf - E^^, E^ - K^,i(Eff), A - iF,,,2(E^^))||e. 

From the Lipschitz continuity (with constant L) it follows 

||k/i - k||lcx=(oxy) < A sup \{xi,yk) - {x,y)\ < L{H Y h), 

{x,y)^€lxY 
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and the same estimate also applies to — /i ^||L~(nxY)- Together with a 
stability estimate for the HMM approximation this yields 

||(Eff - Eh - Kh, 2 (EH))\\e < C{H + /i)||f ||L2(n). 

Splitting the total error E — Eh into the contributions E — Eh and Eh — E 
and using the two estimates, we obtain the assertion. 

Proof of Theorem \5.^A As the terms in the Helmholtz decomposition are orthog¬ 
onal w.r.t. the L^-scalar product, we have 

||Vd||i 2 < IIV0IU2 + ||z||i 2 = ||E - EhWl-. (6.3) 

To estimate z, let (w,wi,ri;2) be the solution of the dual problem 

(w,Wi,'U;2)) = (z,V’)l2 V(?A,' 0 i,'!/' 2), 

and (w//, Wft,, rc/j) the solution of the corresponding discrete dual problem 

Eh{{ipH: V’/i), (wff, vfh, Wh)) = (z, iI}h)l^ i>h, i’h)- 

The analytical and discrete spaces are the same as in the problems for E and 
Eh and therefore not given again here. Because of (V0, z) = 0 it holds ||z|||2 = 
(z, E — Eh)- Thus, it follows 

hWh = (z>E - Eh)l^ = S((eo,61,62), (w,'wi,'u;2)). 

Using the definition of E as exact solution and of Eh as the HMM-approx- 
imation, we deduce 

I|z|Il2 = ^((60,61,62), {w,wi,W 2 )) - H((E,Ki, A:2), {wH,^h,Wh)) 

+ EhiiEn, K?i^i(E//), Kh,2{EH)), (wij, -w/j, Wh)) 

= S((eo, 61,62), (w - vfH, wi - vfh, W2 - Wh)) 

+ {Bh — B){{Eh, Kh,i(EH),Kh,2{EH)), (w/f, Wh^Wh)) 

< C'||(6o,6i,e2)||e||(’W - ’W_f/,Wi -Vfh,W2 - Wh)\\e 

+ C max{||^“^ — ^llL“(axy)) ll«^ — ^^h\\L’^(nxY)] 

II (Ei/, Kh^liEn), Kh,2{E,H))\\e\\i'^H, Wh)\\e- 

According to Theorem 15.21 it holds 

{{{■W -VfH,'Wi - Wh,W2 -Wh)\\e < -f /l)||z||i 2 . 

Hence, together with stability estimates for Eh and wh it follows 

||z||i2 < C{H + /l)||z||i2||(eo,6i,62)||e -f C'?7appr-oa;||z||L2||f||i2. 

To estimate 9, we pose another dual problem: Find w € Hq{V,), W2 S 
H^q{Y)) such that 

Aii^, ^2), {w, W2)) = {9, ^)l 2 V(^, ^2) e X hIo{Y)) 
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with 


-4((V','!/’2), {u,U2)) 


Q JY 


K{x,y){Wip{x) + Vyip2{x,y)) ■ {Vu*{x) + Vyul{x,y)) dydx. 


Again, let us denote by {wH,Wh) S x L‘^{n;W^(Y)) the solution of 

the corresponding discrete dual problem. This dual problem is related to our 
original problem by the equation 

" 02 ), (w, W2)) = -01, V'2), (Vw, 0, W2)) 

for all (^, -!/j 2 ) € x Hj_p(y)) x L^(n; Inserting '0 = 0 

and 02 = 62 , we then obtain 

||0||i2 = A((0,62), (w, W2)) = B{{V 9 , ei, 62), (Vw, 0, W2)) 

= 6((eo, 61,62), {Vw, 0,1x2)) - S((z,0,0), (Vw,0,u;2)). 

By the properties of the Helmholtz decomposition we have {z,'S/wh)l'^ = 0. 
With the same computations as for the dual problem with z, we then derive 

=ve 

||0||^2 = S((eo - 2,61,62), (V(ui - wh), 0 ,W 2 - Wh)) 
-B{{z, 0 , 0 ),{Vw, 0 ,W 2 )) 

+ {Bh - B){^h, 'KhA^H),Kh,2{^H)), (Vrcff, 0, wt)) 

< \A{{9,e2){w- wh,W2 -h)?i))| + C'||z||i21|Vtc + Vy'U)2||L2(nxv) 

+ C max{||/x“^ — ||«: — K^llLoo(nxy)} 

\\{)EH,'Kh,l(EH),Kh,2(EH))\\e\\'^WH + '^yWh\\L^{nxY) 

< C \\\/9 + Vy62||L2(Qxv)l|V(u; — Wh) + Vy(w2 — Wh)\\L^{nxY) 

+ Criapprox\\9\\LA\^\\L'^ + C'll ^ || L2 || 01| L2 , 

where in the last inequality we used the stability estimate for Eh and a stability 
estimate for the solution of elliptic two-scale equations. From a priori error 
estimates for elliptic two-scale problems (see ESI ED]), we know that 

\\V{w - Wh) + Vy(n)2 - w?i)|| < C{H + h)\\ 9 \\L 2 . 

Inserting this, the estimate for the Helmholtz decomposition (lO) . and the 
estimate for z from above, we finally obtain 

||0|li2 < C{H + h)\\ 9 h 2 (||eo|U2(o) + ||V,e2|U2) 


<||(eo,ei,e 2 )|U 

C{H + /l)||(eo, 6 i, e2)||e||^||L2 -|- Cr]approx\\9\\ 


L 2 | 




The estimates for z and 0 together with the a priori error estimate of Theo¬ 
rem 15.21 give the claim. □ 
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6.2 Proofs of the a posteriori estimates 

For the a posteriori estimates we no longer assume higher regularity and there¬ 
fore, need other interpolation operators. 

Lemma 6.4 (Error estimates for the Clement interpolation operator). Let us 
denote by In ■ L^iY) Wj^{Y) the Clement interpolation operator (see fllf) 
with appropriate adaptations to periodic boundary conditions and zero average. 
We define Ih ■ L^{Q] Lfi(Y)) Lfiifil] W( (V)) as 

Ihu{x, y) := 7 h{u{x, ■))iy) 7 x G tl. 

Then the following estimates hold for all u G H^^fY)): 

||m — t/lU||L 2 (TjXSi) < Chi\\V 

||w —xFifc) — \\S/ 

Again, the Clement operator can be defined component-wise for vector functions 
and then fulfills the same estimates. 

Proof. A proof can be found in m- □ 

Lemma 6.5 (Schoberl interpolation operator). There exists an operator Ih '■ 
Ho(curl) —> VjjQ such that for every u G Ho(curl) there exist 0 G Hq{LI) and 
z G nun) with 

u — IhVl = V9 -I- z. 

The decomposition fulfills the estimates 

H~^\\d\\L‘^{Tj) + l|V0|h2(r^.) < C\\u\\L2(^cjTfi^ 

+ l|Vz|h2(T.) < C||curlu|h2(,i^.). 

Together with the trace ineguality we moreover have the estimates 

||z||l2(f,.,) < Ci7j/^||curlu|h2(^.^.), 
where, of course, Tj can also be substituted by Ti. 

Proof. For the construction of Ih and a proof of the estimates we refer to [5S] . 
Additional details on Ih can also be found in [531133] ■ CH 

Proof of Theorem \ 5 . 5 [ First of all, we derive an error identity. From the defini¬ 
tion of the error terms and Proposition 14.31 we deduce 

^((60,61,62), {if, if 1,1(2)) 

= S((E,Ki,iF2),(^,tAi,V’2)) 

- Bi(EH, KhpC^H), Kh,2i'EH)), (ip, ^/’i, V’2)) 

= (f, if) - B{{Eh, K^,i(Eff), iF^,2(Eff)), {if, iP„iP2)) 

- Bh{{^H,n.hp(EH), Kh.A'En)), (tp - tpH^ V’l - '02 - i’h)) 

+ Bh{{^EH,KhA'EH),Kh^ 2 {'EH)),i^P,'^Pl, 1 p 2 )) 

- Bh{{{'EH,Kh,i{^H), KhafEn)), {ipn^ i’h)) ■ 

' -V-' 
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With the definition of the residual (lO) . this gives the following error iden¬ 
tity for all e Ho(curl) x x {n-, HIq{Y)) and all 

e X L\^-Wl{Y)) X L^{n-Wl{Y))-. 


S((eo, 61,62), (■»/), ■0i,V’2)) 

= -{ReshCEn, Kh^iiEn), Kh, 2 {EH)), (xp - xpfj, ' 0 i - ' 0 / 1,02 - xph)) ( 6 . 4 ) 
+ [Bh - B){{Eh, K„,i(Ei^), Kh,2{EH)), (0,01,02)). 

We choose 0 = cq, 0 i = ei, 02 = 62 and xpfj = Ineo, xpf^ = hei, xph = Ihe2 
in the error identity ( 16 . 41 ) with the interpolation operators Ih from Lemma 16.51 
and Ih from Lemma [ 6.41 Using the coercivity of B, we obtain 

||(eo, 61, e2)||e 

<C(|(Res„(Eif,K„,i(E^^ ), R^/i.2(E_ff)), (60 — Ih^o, 61 — //iCi, 62 — 4/162))! 

'-V-' 

:=I 

+ \{Bh- B){{Eh, KhA^H),Kh,2{^H)), (eo, 61,62))!). 

'-V-' 

:=II 


To estimate I, we insert the decomposition of eo — Ih^o = V0 + z according 
to Lemma 16.51 and thus obtain 


1= ^^^(curlE//-I-curlj^K/ip) • (curlz*-I-curly(ei -//iCi)*) 

JQ Jy 

-I- divy K/i_i divy(ei - luei)* 

- Kh{EH + yyKh, 2 ) ■ {Ve* +z*+ Vj,(e 2 - Ihe 2 )*)dydx 

- [ i{V9*+z*)dx . 

Jn 

Integrating by parts locally, inserting in — fn and noting div f = 0 yields 
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divy{Kh{EH + ^yKh,2)) ■ (62 - Ihe2)* dydx 


jeJ -'j iei 


-EEL {khCEh + VyKh, 2 )) • n (e2 - 462)* dadx 


j&J iel 


IdS; 


-E E/ {{h + i - in + khCEh + VyKh,2)) ■ z* dydx 


j&J iel ' 


As Kh are constant on the cells Tj x Si, the correctors Kh are constant with 
respect to x and linear with respect to y, and 'Em is linear with respect to x, 
all terms with two derivatives with respect to the same variable cancel out. We 
derive by a rearrangement of sums and the Holder inequality 


• z* da 


1 = 

E / 

/ (/r^ycurlE/f + curly K?,_i)) X ndy 



L7v 


- S L,.. 


e*da 


{Kh{EH + VyKh,2)) -ndy 
^ LL ^(curlE/f + curL K/j^)) x n + divj^ 

3eJ{i,k)GS(Thy'^^ 

• (ei — IhS-i)* dadx 

^ / /- [{^^h{EH y'^yKha)) ■VL]p.^{e2-Ihe2)* dadx 

-f / drvxi I Kh{EH+ '^yKh, 2 )dy\ 9 * dx 
/ (i — iH+tH+f HhiEn + '^yKh,2) 




E 

{hl)^£{TH) 


^JY 


z* dx 


(/r^LcurlEij+ curlj,K;,4)) xndy ||z|L2(^^,) 




^JY 


{Kh(EH + VyKh,2)) -ndy \\ 0 \\LyF,ii 

J Fji 


^Fik 


+ E 

U,i)&£irH) 

+ E E (/x^ycurlE// + curly Kh,i)) x n + divyK/^^n 

jsJ {i,k)(iS(Th) 

||(ei — -^/tei)|L2(jvxAfc) 

-f E! E! ||(e 2 — 

jsJ {i,k)(iS(Th) 

+ E||'^^^^(/ + VyA:/i,2)dy) II^'IIl2(t,) 

+ E(|I^ “ y dy ) ||z||L2(Tj)- 


Using the estimates for the interpolation operators from Lemmas 16.41 and 16.5 
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we get 


I < 


3 1 


U,1)&S(Th) 


^JY 


+ curly Kh,i)) xndy 


Fji 


I curleo||L2(^ 


+ E I f / + ^yKha)) -ndy] ||eo|U 2 (^ ) 

III,, j 


{3,1)&S{Th) 


'-JY 


+ E E I (curl E// + curly K?,,i))xn + divj^K,,,in]^.^ 

3eJii,k)eS{Th) 

W'^yClWL'^iTjXujp^J 

+ E E Chli^inniEH + VyKha)) ■ n]f,J| \\Vve2\W(T,xui,j 

jeJ {i,k)eS{Th) 

+ Y, CHj div^^/ Kh{^H+ ^yKh,2)dy'^ \\eo\\L'^{a,T^) 

+ ^OTy(||f-fff|| + ||fo+ / Kh{^H + VvKh,2)dy\^ 

jeJ 

II curleo||L2(^ ). 


Applying the Cauchy-Schwarz inequality, we obtain the desired local estimators 
and terms like X^jgJ II curled| l 2 (£;;,^ ). As the triangulation is shape regular, 
each element Tj only appears in a finite number of these neighborhoods and 
this number can be bounded above by a uniform constant (independent of H, h). 
Thus, we derive 


I<C'( E ’lAi) ^ ||curleo||L2 (o)+c( Y ^ l|eo||L2(o) 

{3,i)(i£(TH) U,i)&£iTH) 

1 /2 

+c(i: i: l|Vyei||L 2 (nxv) 

jGJ ii,k)e£{Th) 

1^2 

+c(i: i: ^i,ifc,2) l|Vye2||L2(nxv) 

jeJ ii,k)e£{Th) 

+ ^'{j2^h) l|eollL=(n)+t^(E Cj+4i) ^ l|curleo||L 2 (n). 

3&J 3&J 


Of course, all norms of the errors eg, ei, and 62 can simply be estimated by 
||(eo, ei, e 2 )||e- 

To estimate II, we just split the integral into local terms and use the Holder 
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inequality 


11 = 



Vi) - H ^(x,y))(curlEff+ curlj^K,i_i) 


• (curl Cq + curlj, e*) 

+ - k{x, y))(E// + VyKh^) ■ (gq + '^yS-Vl dydx 

< ll(eo,ei,e 2 )||e. 

iGJ ie/ 

Dividing each estimate by ||(eo, ei, e 2 )||e and combining both gives us the a 
posteriori error estimate. □ 


For the proof of the lower bound we need local bubble functions. Let us 
denote by At,;, / = 1,..., 4, the barycentric coordinates of a tetrahedron T and 
by At,;, I = 1,2,3, the barycentric coordinates of a face F. The local bubble 
functions on elements and faces are defined as 

4 3 

^Jt := 256]^ At,;, := 27]^ At,;. 

;=i 1=1 


They fulfill 0 < < 1, suppV'T C T, and supp^'T C u)f- We also define 

a continuation operator Pp : L°°{F) —>• L°°{u}p) as the constant extension of a 
function in the direction perpendicular to the face F, for details see [38]. The 
following inequalities can be proven with standard scaling arguments and the 
properties of the bubble functions (see [351 Proposition 3.37] for details and the 
proof). 

Lemma 6.6 (Inverse inequalities). For all g and all tetrahedra T it holds 

Mh^T) < C\{g,ipTg)T\, 

IIV^Tff||L 2 (T) < C'||g||T 2 (T), 
l|V(V’Tff)||T 2 (T) <Cdiam(T) ^|| 5 ||t 2 (t). 

Furthermore, for all f gF^\f and faces F it holds 

ll/lli2(T) <C|(5 ,V'fPf(/))f|, 

IIV'fPf(/)||t^(..) < Cdiam(E)i/2||/||T2(T), 
I|V(V'fPf(/))||t 2 (^,) < Cdiam(E)-i/ 2 ||/|U.(^). 

Proof of Theorem [13 First of all, corresponding to the error terms we intro¬ 
duce the following functions: 


^i,i{x) =+ J Kh{-,y){'EH + '^yKh^2i-,y))dy^ix), 

Wj;,i(a;) = V'Ft(2;)Pft ( / (/r^^^curlE//-P curly K/i,i)) x ndy )(x), 

JFji/ 


'j,ik,i{x,y) = XTj(a:)V'F,-, (y) 


-^F (curlE f + curly K,it)) x n-tdiVyK,,,in]^ )(?/), 
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Wj,2{x) = tpTj{x)divj;(^J Kh{-,y){'EH + S/yKh,2{-,y))dy^{x), 

Wji^ 2 {x) ='ipF^i{x)PFji( / Kh{-,y)(EH+ VyKh, 2 {-,y)) ■ ndy )(a;), 

^Fji/ 

WJ,^k, 2 {x,y) = XTAx)'iljp.^{y)Pp.^{[Kh(EH+ VyKh, 2 ) ■ n]p,J{y), 

where XA denotes the characteristic function of the set A. 

Now the error indicators can be estimated separately: By partial integration 
they can be interpreted as the residual tested with the localized functions. The 
error identity and Lemma 16.61 are then used each time to bound the indicators 
by the total error. 

As Eff is linear and Kh,i constant with respect to x, we get with integration 
by parts and suppWj_i C Tj 

0 = / / curl3;(/x;;;^(curlE//+ curlj^ K,i_i)) • dyda: 

dn JY 

= / / (curl Ejj + curly K/j^i) • curl w* dyda;. 

JnJY 

Therefore, we obtain for the first error indicator 


!4i<c 


J (fn + j + 




yJAh,2 


dx 


= C | —(Res?i(E//, Khx,Kh,2), 0,0)) + {in — f, 'Wjx)L^ \ ■ 

If we choose ijj = Wj^i, = 0, '!/'2 = 0 in the error identity (16.41) . we get with 
Lemma 16.61 

I - {Reah{EH,Khx,Kh,2), (wj,i,0,0)) + {in - 

= |i 3 ((eo,61,62), (wj,i,0,0)) + (6 - Bh){{EH (wj.i,0,0)) 

+ (fo - f,Wj,l)i 2 | 

< (^11(60, 61 , 62 )||e(T 3 x>")ll’*^J.lllH(curl.Ti) + C'llfe “ f II L2(t, ) || Wj.i || i2( jv ) 

1 y/2 

+ 11^1.1 llH(curl.ri) 

iei 

< C\\iH - t\\L^{Tj)H~^ r]j^i + CTL^^IKeo, 61, e 2 )||e(Tixi') ? 7 i,i 
1/2 


+ ch-^(J2c 


iei 


Vj,i- 


All in all, after multiplying by iL| rjj l, this gives the local estimate for 
For r]j ^2 we get with the properties of Wj ^2 and an integration by parts 




< C 


w* 2 dx 


J divj;(^J Kh(EH + VyKh,2 
-J (^J Kh{EH+ '^yKh,2)dy^ ■Vw*2dx 


= C 

iTj '^JY 

= C \{Resh{EH, Khx,Kh,2), (Vwy_2,0,0))|, 
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where the last equality holds because of curlV = 0 and {i,Vwj^ 2 )L^ = 0- Now 
we can use the error identity (16.41) with tp = S/wj^ 2 , '0i = 0, -02 = 0 and obtain 
with Lemma 16.61 


|(Res;j(E//, Kh,i,Kh,2), (Vri;j,2,0,0))| 

= I - S((eo, ei, 62), (Vzc,-2,0,0)) + {Bh - S)((Eff, Kh,i, Kh,2), (Vu;,,2,0,0))| 

1^2 

< C'(^||(eo,ei,e2)||e(T3x>")ll^'^i.2||L2(T,-) + 

iGl 

1/2 

- ^ll(eo, ei, e 2 )||e(Tjxy) ^j .2 + Hj ^(X‘’fi) ■ 

i&I 


Multiplication by iL| rj- ^ gives the local estimate for r]j^ 2 - 
For rjji^i we have the estimate 


V]i,i 

Hu 


< c 

/ 

/ (/r.^(curlE//+ curLK,j 1)) X ndy 

-wilder 


Jfu 

LJy 

Fji 


An integration by parts and the linearity of E/f yield 


(curlE// + curlj^K,i_i)) x ndy ■vf*i^da 
iFu^JY ^Fji 


= / ^(curl E//+ curlj, K/i_i) • curl dyda: 

diip., Jy 

= (Res/t(Ei/,K?i,i, A:;j, 2), (wj74,0,0)) + [ (f - fo) • wh ^ 

+ J + J khCEh+ yyKh,2)dyj -Wji^^dx. 


dx 


The integrals over can be split into two element integrals over Tj and T;. 
In the third term we recognize the known error indicator 1 and in the second 
term the data approximation error indicator Q. Using first the error identity 
(16.411 with tp = xpi = 0, and ip 2 = 0 and Lemma 16.61 we obtain 


Hji 


— II ) ^ 2 ) II e(a;F 


1/2 


V)||Wj7,l||H(curl.LDj=.^.,) 
Ilwj-jp ||H(curl,iiF^.j) 


+ l|f - l|L2(Tj)||Wji,l||L2(7^) + ||f - {H\\L^{Ti)\\YVjl^i\\L^(Ti)^ 

1/2 

< c(^\\{eo,ei,e2)\\e{u,p.^xY)H~i^ + H~l^(^Cj^ + Cfi) Vji.i 


+ {H- ^ rjj^i + Hi ^ 77 / 4 ) + {H- ^ Q + Hi ^ 0)Vji,i^ ■ 


Due to the regularity of the triangulation the quotients Hji/Hj and Hji/Hi can 
be bounded above and below. Thus, multiplication by Hjir]~i\ together with 
the already derived estimate for rjj^i yields the desired estimate for ryjzp. 
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For r]ji ^2 we obtain by the properties of Wji^ 2 , an integration by parts, and 

(f, VWj7,2)l2 = 0 




H 




< c 

[ 

/ {Kh{EH + yyKh,2)) ■ ndy 

w*i2da 


dFil 

Uy 



= C 


= c 


J khCEh+ yyKh,2)dy^w*i2dx 

f / / khCEh -\-VyKh,2) ■'^'W*i2dydx 

jcjf-, 'jy 

J divjj^y Kh(EH + yyKh,2.) dy^w*i 2 dx 

<^Fji y 

— (Resh^En, Kh,2), (Vwji,2,0,0)) 


When split into integrals over Tj and Ti, the first term can be identified with 
r]j ^2 and rii^ 2 , respectively. For the second term we can insert again the error 
identity (16.41) with ■0 = Vwji, 2 , 0i = 0 = 02- With Lemmawe then get 


^1,2 


< C'(||(eo,ei,e2)||e(<iF^.jXV)||Vw;j7,2||L2(tDF^.,) +Hj Nj.2||'Wj7,2||l2(Tj; 

1^2 

+ Hi ^ Vl,l \\'>J^jl,2\\L'^iTi) + Cji + Cfi) I|VWj7,2||l2(cDf^,)) 


< C'dKeo, ei,e2)\\eiu;p., xY)Hji^ r]ji^2 + {H^ ^ 77^,2 + H^ ^ r]i^2) Vji,2 

i 

Together with the already derived estimate for 77 ^^ this gives us the local esti¬ 
mate for 7777 , 2 ■ 

For 777,i7c.i we have the estimate 


hik 


< c 


'Ti J Fit 


[(/r^^(curlE_f/ -|-curl^ K/j,!)) x n-|-div^ K 7 ,,in]p.. 


• dadx 


With an integration by parts and the linearity of with respect to y we 
obtain 

/ / [(M0^(curlE//+curlyKh,i)) X n-|-divyK;i,in]7^,.^ ■ ;^dcrda; 

dTj JFik 

= /^^^(curlE//-I-curly K7i,i)-curlj^ w*^_;^-l-divj^ K/i,i divj^ dj/da; 

JTj ddlp_^ 

= {ReSh(E,Hi Kh,2), (0, W7,ife,i, 0)). 

Inserting the error identity (j6.4|) with 0 = 0 = 02, 0i = 'Wj,ik,i and using 
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Lemma l6.6l we obtain 



For r]j^ii ~^2 we derive with integration by parts and the linearity of Kh ^2 with 
respect to y 



cl —(Res;i(Eif, K/i4, Kh,2), (0,0, Wj^ik,2)) 


Inserting once more the error identity this time with ip = 0, ipi = 0, and 

^2 = Wj,ik, 2 , we obtain with Lemma l6^ 



This gives us the local estimate for The global estimate now follows by 


summing up the local estimates. 


□ 


Conclusion 

In this paper, we suggested a new Heterogeneous Multiscale Method (HMM) 
for the time-harmonic Maxwell equations. The basis is a homogenization re¬ 
sult for a curl-curl-problem obtained with two-scale convergence. Divergence- 
regularization is applied to the corrector for the curl and thus we can get rid of 
the divergence-free constraint. The HMM can bee seen as direct finite element 
discretization with numerical quadrature of the two-scale homogenized equa¬ 
tion, which is the crucial observation for the numerical analysis. The a priori 
error analysis shows that the HMM converges with the same rates as finite ele¬ 
ment methods for curl-curl-problems without oscillating coefficients (Theorems 
15.2115.3p . The a posteriori error estimators are reliable and efficient (Theorems 
[531 Ol) and can be used for adaptive algorithms future work. 
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